The sky brightness and transparency in z-band at Dome A, 

Antarctica 

Hu Zou 1 ' 2,3 , Xu Zhou 1 ' 3 , Zhaoji Jiang 1 ' 3 , M. C. B. Ashley 4 , Xiangqun Cui 3 ' 5 , Longlong 
Feng 3,6 , Xuefei Gong 3,5 , Jingyao Hu 1,3 , C. A. Kulesa 7 , J. S. Lawrence 8 ' 9 , Genrong Liu 5 , 
D. M. Luong-Van 4 , Jun Ma 1 , A. M. Moore 10 , C. R. Pennypacker 11 , Weijia Qin 12 , Zhaohui 
Shang 13 , J. W. V. Storey 4 , Bo Sun, 12 , T. Travouillon 11 , C. K. Walker 7 , Jiali Wang 1 ' 3 , Lifan 
Wang 3 ' 6 , Jianghua Wu 1 , Zhenyu Wu 1 , Lirong Xia 5 , Jun Yan 1 ' 3 , Ji Yang 6 , Huigen Yang 12 , 
Yongqiang Yao 1 , Xiangyan Yuan 3 ' 5 , D. G. York 14 , Zhanhai Zhang 12 , Zhenxi Zhu 3 ' 6 

ABSTRACT 

The i-band observing conditions at Dome A on the Antarctic plateau have 
been investigated using data acquired during 2008 with the Chinese Small Tele- 
scope ARray. The sky brightness, variations in atmospheric transparency, cloud 
cover, and the presence of aurorae are obtained from these images. The median 
sky brightness of moonless clear nights is 20.5 mag arcsec -2 in the SDSS % band 
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at the South Celestial Pole (which includes a contribution of about 0.06 mag 
from diffuse Galactic light). The median over all Moon phases in the Antarctic 
winter is about 19.8 mag arcsec -2 . There were no thick clouds in 2008. We model 
contributions of the Sun and the Moon to the sky background to obtain the re- 
lationship between the sky brightness and transparency. Aurorae are identified 
by comparing the observed sky brightness to the sky brightness expected from 
this model. About 2% of the images are affected by relatively strong aurorae. 

Subject headings: methods: data analysis — methods: statistical — techniques: 
photometric — telescopes 



INTRODUCTION 



In selecting observatory sites for ground-based optical/IR astronomy, some of the most 
important considerations are the night-sky brightness, transparency, seeing, number of clear 
nights, humidity, and photometric stability. The Antarctic plateau offers some attractive 
advantages for ground-based astronomical observations. Site testing over the past decade 
has revealed that Antarctica, relative to temperate latitude observatories, has lower in- 
frared sky brightness, better free-atmosphere seeing, greater transparency, a l ower turbulent 
boundary layer, and much lower water vapor conte nt (see, e.g., reviews by lAristidi et al. 



f l2005al ) JBurton et all (12005( 1. and lStorev et al.l (120071 ) ). About half the area of the Antarctic 
continent is at an elevation of over 3000 m above the sea level, and the year-round average 
temperature on the plateau is about — 50°C, dropping to below — 80°C on occasion. The 
aver age of the water va por column over the plateau has been predicted to be less than 250 
fim (IChamberlinl 120011) and h as been measured at Dome A to be below 100 /xm for 25% 
of the time (jYang et al.l l2010l ). Wind speeds at the summits of the plateau are very low 
close to the surface and the median thickness of the turbulent boundary layer has been mea- 
sured at Dome A as 13.9 m (jBonner et al.ll2010l ). These exceptionally favorable geographical 
and meteorological conditions result in substantial advantages for millimeter, sub-millimeter, 
infrared, and optical astronomy. 

Several large astronomical facilitie s are already oper ating at some sites on the Antarctic 
plateau, e.g., the South Pole Telescope (IRuhl et al.ll2004l ) at the US Amundsen-Scott station. 
Proposed facilities inc lude a 2.5 m optic a l telescope, PILOT (n ow called PLT) at Concor- 
dia Station, Dome C (jStorey et al.l 120081 : lLawrence et al . 2009a ), and a 2 m class infrared 
telescope planned for deployment to Dome Fuji by Japan (ITakato et al.l 120081 ) . 



At Dome C (3250 m elevation), the average seeing is 0."27 above a ~ 30 m ground 
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layer, and the average surface wind speed is 2.9 m s _1 ( [Lawrence et al.l 12004 ; lAristidi et al. 
2005aU bl). Dome A is the highest point of the continent, with an elevation of 4093 m. It 
might be reasonably predicted that Dome A could be as good as or even a better astronom- 
ical site than Dome C , with better seeing, higher transparency, and thinner surface layer. 
Saunders et al.l (120091 ) compared the sites Dome A, B, C, F and Ridge A and B in their cloud 
cover, free-atmosphere seeing, precipitable water vapor, temperature, and auroral emission, 
and concluded that, overall, Dome A might be the best of the existing bases for astronomical 
observations. 

Following the success of the first Chinese expedition team to Dome A in 2005 January, 
the Chinese Center for Antarctic Astronomy (CCAA) began planning for the dep l oymen t 
of a telescope array called Chinese Small Telescope ARray (CSTAR; lYuan et al.l (120081 )) 
described in the following section. 

In this paper, we use CSTAR images obtained during 2008 to assess the night sky 
brightness at Dome A. Both the variation in atmospheric transparency and cloud cover are 
also derived. We then analyze the contribution to the sky brightness from the Sun and the 
Moon and investigate the correlation between sky brightness and atmospheric transparency. 
Finally, the contribution from aurorae is assessed. 



2. INSTRUMENT AND OBSERVATIONS 

The CSTAR telescope array was built by the Nanjing Institute of Astronomical Optics & 
Technology (NIAOT) and consists of four co-aligned Schmidt telescopes on a fixed mounting 
pointing toward the South Celestial Pole. Each telescope has a different filter: g, r, i, and 
open. The detectors are Andor DV435 lkxlk CCDs with 13 /im pixels, giving a plate 
scale of 15" pixel -1 . The entrance pupil diameter is 145mm (effective aperture: 100mm) 
and the field of view is about 20 deg 2 . The scientific goals for CSTAR were primarily 
the detection of variable stars and measurements of sky brightness, the stability of the 
atmospheric transparency, and cloud cover. 

The CSTAR telescopes, cameras, and computer system were tested under conditions 
of very low temperature, low air pressure, and long-duration continuous operation. The 
observational testing of the arra y was done a t theXinglong Station of National Astronomical 



Observatories of China (NAOC; IZhou et al.l (J2010aJ))- 



In 2008 January, CSTAR was transported to Dome A and successfully installed, with 
the intention of making observations of the South Celestial Pole region from 2008 April to 
August. Power, heat, and Internet connectivity for CSTAR were provided by PLATO (the 
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PLATeau Observatory ), developed by the University of New South Wales ( lYang et al.l 12009 



Lawrence et al.ll2009al ). 



After deployment to Dome A, only the z-band telescope (CSTAR^l) operated well, ob- 
serving throughout the winter with just a few short interruptions; the other three telescopes 
had either small problems with their optical alignment or their control computers. How- 
ever, all four telescopes operated sufficiently well to allow us to derive some information on 
sky conditions in the other filter bands, and this analysis will be published in a subsequent 
paper. Over 310,000 images were obtained, with integration times of 20 s or 30 s, between 

200 8 March 4 and 20 08 August 8, giving a total exposure time on the sky of about 1728 
hr (IZhou et al.ll2010bl ). While CSTAR was observing, roughly one-third of the catalogs and 
one original image were transferred via Iridium communication satellite each day. All the 
catalogs and images were brought back by the 25th Chinese Research Expedition Team, in 

2009 April. The catalogs contain the ambient temperature, observing date and time, expo- 
sure time and background brightness, plus the instrumental magnitudes in three different 
apertures for all detected point sources calculated by an automatic aperture photometry 
pipeline. For each image about 10,000 sources down to % ~16 mag were detected. Several 
variable stars were identified and the complete observational catalog has now been released 
(IZhou et al.l l2010bl ). As the catalog has continuous observations over several months, the 
data are ideal for the detection and monitoring of variable stars (J. H. Wu et al. 2010, in 
preparation). 

Observations began on 2008 March 4 and ended on 2008 August 8, a few weeks earlier 
than expected, due to a power failure. Images were obtained whenever the Sun elevation was 
sufficiently below the horizon. There were two time periods during which no observations 
were made due to issues with the computer system: 2008 May 19-29 and 2008 June 30 to 
2008 July 15. 



3. RESULTS AND DISCUSSIONS 



3.1. Sources of Sky Brightness 



Night sky brightness is one of the key parameters for qualifying a site for ground-based 
astronomical observation. Many sources contribute to the sky background including scatter- 
ing from the Sun and the Moon, airglow, zodiacal light, aurorae, star l ight and interstellar 
dust scattering, extra galactic light, and artificial light contamination ( jLeinert et al.l 119971 ; 
Benn fc Ellison! 1 1998al ) . The contributions of th ese sources to the sky brightness at Dome C 
have been reviewed by iKenyon Storey! (120061 ) and this is generally applicable to Dome A. 



-5 - 



Artificial light pollution is essentially non-existent in Antarctica. The main contribution 
to the sky b ackground is usually fro m the atmospheric scattering of the light of the Sun and 



the Moon. iBenn fc Ellison! ( I1998af ) showed that zodiacal light (i.e., sunlight scattered by 
interplanetary dust) may contribute up to half the intrinsic sky brightness when the Sun 
and the Moon are down. However, this is not an issue for CSTAR observations since the 
field that we are observing is centered on the South Celestial Pole, with ecliptic longitude 
A of 270° and latitude (3 — 66. °6, and is far from the region affected by zodiacal light. On 
the other hand, we expect an increased contribution from the Milky Way, since the Galactic 
latitude b of the South Celestial Pole is —27°, relatively close to the Galactic plane. We 
quantify this effect in Section 13.31 

Another important source of light is from the aurorae, which are common in polar re- 
gions. Aurorae are generated as electrons from the solar wind which are trapped in the 
upper atmosphere by the magnetic field of the Earth. The intensity and frequency of oc- 
currence of aurorae depend upon the solar activity. From satellite images and ground-based 
observations, aurorae are known to be con centrated into an oval -shaped annulus centered on 



the geomagnetic south pole. As noted by lDempsey et al.l (120051 ). the more intense the solar 



activity, the larger the extent of the oval. Domes A and C both lie within the "hole" in the 
center of this oval-shaped region and so are less seriously affected by aurorae than sites lying 
on the oval itself. In 2008, solar activity was close to its minimum level so that, perhaps 
paradoxically, Dome A would have experienced more aurorae than in an average year. The 



maximum of solar cycle 24 is predicted to appear in 2012 December (IKilcik et al.ll2009l ). so 



the number and intensity of aurorae appearing above Dome A should then be at a minimum. 

To better understand the effect of the Sun, the Moon, and cloud on the sky brightness, 
we modeled the contributions of the Sun and the Moon light without the influence of clouds 
in Section 13.61 and 13.81 After removing these two sources of sky light, we consider the 
relationship between the residual sky background and atmospheric transparency (see Section 
13. 9p . The remaining variable component of the sky brightness can then be attributed to 
aurorae (see Section [3. 10p . 



3.2. CSTAR Measurements of Sky Brightness 

The catalogs processed by the CSTAR pipeline contain the sky background values in 
analog-to-digital converter counts (ADU). The instrumental magnitude, m, in the CSTAR 
photometric system is defined as 



mcsTAR = -2.51og(ADU) + 25. 
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Zhou et al.l (l2010bl ) used the USNO-B1 (U.S. Naval Observatory) star catalog to flux cali- 
brate the CSTAR images and derived a calibration constant o f icsTAR — jnsNQ = 4.16±0.12. 



SDSS i-band magnitudes are very close to AB magnitudes (IBohlin et al.l 120011 ). so we can 
transform o ur z-band magnitud es into a top-of-the-atmosphere flux. The magnitudes in the 
AB system (jOke fc Gunnlll983l ) are defined as follows: 



-2.51ogF„ -48.6, 



where m is the AB magnitude and F v is the flux of the specified wavelength band in the 
units of erg s -1 cm -2 Hz -1 . Then we can convert from CSTAR raw ADU counts (normalized 
to a 20s exposure) to an z-band sky flux per square arcsec using 

F = ADU x io- 27 798 /15 2 = 0.708 x 10 -30 ADU, 

where we have assumed that the typical atmospheric extinction in clear photometric nights is 
about 0.05 mag/airmass at i, taken from the observational results of the Carlsberg Meridian 
Telescope (CMT) at La Palma, and 15 2 is the pixel area in square arcseconds. 

The surface fluxes in the rest of this paper are expressed in units of 10 -30 ergs s -1 cm -2 Hz -1 arcsec 

Figure [I] shows the sky brightness variations throughout the whole observing period of 
2008. A diurnal cycle is evident even in midwinter. At such far-southerly latitudes as Dome 
A, the Moon is always fairly full whenever it is above the horizon from April to August, 
creating a strong correlation between lunar elevation (including the correction for parallax) 
and sky brightness. 

Figure [2] shows histograms of the sky brightness distribution at Dome A. The most 
probable sky brightness, across all lunar phases (left panel of the Figure), is found to be 
20.1 mag arcsec -2 and the median is 19.8 mag arcsec -2 . About 80% of the images have sky 
backgrounds darker than 19 mag arcsec -2 . 

The right panel of Figure [2] shows the sky brightness distribution of images taken only 
on moonless and clear nights during 2008 June. In these images, the Sun and the Moon 
elevation is below —18° and there is no cloud. We find that the median sky brightness is 
20.5 mag arcsec -2 . When it is dark but cloudy, the sky brightness increases to about 20.2 
mag arcsec -2 . When the Moon is full and it is clear, the brightness is around 19.1 mag 
arcsec -2 . When the Moon is full and there are clouds, the average sky background is 17.9 
mag arcsec -2 . 
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Moles et al.l (120101 ) have summarized the night sky brightness in UBVRI measured at 



vario us observatories during dark time. We used the transformation equations (IJordi et al. 



20061 ) between SDSS magnitudes and UBVRI to compare our i-band measurements with 



-7- 



the sky backgrounds of La Palma, Cerro Tololo, and Paranal with Dome A. These trans- 
formations, which are appropriate for stars, are slightly different from those that should be 
used for the sky background because of the presence of strong emission lines. We calculate 
the followi ng median i-band sky b rightnesses: 20.10 mag arcsec -2 at La Pal ma (at sunsp ot 

'aranal (at s unspot max imum; iPatatl ( 120031 )). 



Benn & Ellison! (I1998M )). 19.93 at 



Walker. A. I (1l987l . ll98Sh ) and 19.57 at Calar 



minimum; 

20.07 at Cerro Tololo (at su nspot minimum ; 
Alto (at sunspot maximum; ISanchez et al.l (120071 ) ) . The CSTAR measurements reported 
here, with a clear dark sky brightness of 20.5 mag arcsec -2 , support the tentative conclusion 
that under moonless clear conditions, Dome A has a darker sky background than the above 
astronomical sites, even allowing for calibration and transformation uncertainties of up to 
several tenths of a magnitude. 



3.3. Correction for Galactic Background 



To compare our observations of the sky background at the South Celestial Pole with 
observations taken well away from the Galactic Plane we need to estimate the contribution 
from diffuse Galactic light. The Pioneer 10 and 11 spacecraft made observations of the total 
Galactic plus extragalactic sky background in blue (395-485nm) and red (590-690nm), from 
beyond the asteroid belt . The results for selected fields, including the South Celestial Pole, 
are given in lToller et al.l ( 119871 ). The Pioneer values have a resolution of about 2° on the sky 
and include the light from stars fainter than about V ~ 8. We extrapola te the Pionee r red 



data into the i band by using the ratio of line-of-sight integrated starlight (jMattiralll980l ) and 
the ratio of line-of-sight extinction at the different wavelengths. This leads to an estimate of 
26 ±8 x 10 -10 erg cm _2 s -1 sr -1 nm -1 , which is a 23.75 mag contribution to our median of 20.5 
mag (i.e. ~ 5% of the total sky background). In other words, we would expect the median 
sky brightness in regions well away from the Galactic plane to be a further ~ 60 mmag 
fainter. 



3.4. Transparency Variations 

Starlight is attenuated as it passes through Earth's atmosphere due to molecular ab- 
sorption and scattering from molecules and aerosols. Ideally, there would be stable, dry 
air above the site and, in the absence of clouds, the observed flux of each non-variable star 
would show little change. Normally, it would be possible to derive the absolute extinction 
for the site by observing the same star at two very different airmasses or by observing stars 
of known brightness that happen to be at two different airmasses. However, in the case of 
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CSTAR, with only a 4.°5 field, neither technique is practical. We thus confine our discussion 
to variations in transparency, with the absolute value of the extinction yet to be determined, 
although presumably it is low at Dome A due to the altitude of the site and the low aerosol 
content of the atmosphere. 

To assess the variations in transparency, we compared the observed brightness of an 
ensemble of bright, unsaturated, stars to the brightness of these stars on a reference image, 
chosen as one containing the greatest number of detected point sources. All images were 
scaled to the same exposure time, 20 s. 

Figure [3] presents the time variation of atmospheric transparency during the periods of 
darkness throughout 2008, while Figure H] shows an expanded version covering 10 day periods 
in each month from April to July. From April 15 to June 2, the transparency was relatively 
stable with many days showing high transparency. In June and July, the sky reaches its 
darkest due to the low elevation of the Sun; the transparency became a little worse in July. 
The median values of the excess extinction, above that of our nominally clear reference 
image, for the months from April to July are 0.056, 0.016, 0.127, and 0.49, respectively. 

Figure [5] shows the distribution of the transparency. Most of the images have only 
a small increase in transparency relative to the reference image: 90% of the images have 
less than 0.7 mag of transparency, 80 % have less than 0.4 mag, and more than half have 



less than 0.1 mag. Ilvezic et al.l ( 120071 ) show that in conditions of up to 1 mag of reduced 
transparency, photometry is limited by photon noise, not by variations in the cloud extinction 
itself. Therefore, conditions were suitable for accurate photometry for > 90% of the time 
during 2008 at Dome A. 



3.5. Percentage Cloud Cover 

Here, we use our observed relative transparency variations to infer the distribution of 
the optical thicknesses of cloud at Dome A during 2008. We note that any reduction in 
transparency could be due to combinations of cloud, the atmosphere, possible ice formation 
on the entrance pupil of the telescopes, and the exhaust plume produced by the diesel 
generators (about 45 m away from CSTAR). Our results therefore represent upper limits to 
the cloud. We expect window frosting to be negligible, due to the design of the front window 
of the telescope which is heated all the time while observing. Table [2] lists the percentages 
of images in various extinction ranges — chosen to correspond to the flux percentage of stars 
relative to the flux in the reference image as shown in the first column. About 51% of the 
images were obtained in excellent conditions with excess extinction from any cause of less 
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than 0.11 mag. 

For a rough comparison with another observatory, Table [3] shows the approximate frac- 
tion of cloud of various thicknesses at Mauna Kea. The data for the table come from the 
Gemini Observator}0 and are based on nightly logs from the United Kingdom Infra-red Tele- 
scope over a 10 year period. The relative transmission levels Gemini provides are in l/-band 
magnitudes (these levels are subjectively assumed by them for the purpo se of integration 



time calculation). We use the calculated atmospheric extinction (jKinglll985l ) to convert these 
V-band extinction in clear weather to i band by simply applying the differential extinction 
between these two passbands (V — i = 0.07 from the observational results of the CMT for 
La Palma on Web site of |http: //www. ast . cam. ac .uk/ ~d we/SRF/camc_ext inction.html ). 
Here, V — i — 0.07 is just like a systematic zero point to transform the cloud extinction 
ranges of different cloud categories between these two bands, although it is unclear about 
how the cloud extinction varies in different passbands and from site to site. Table [3] shows 
that the cloud cover at Dome A falls into the "cloudy" or worse category for a mere 2% of 
the time, compared with 30% at Mauna Kea. There are patchy clouds at Dome A for about 
31% of the total usable time, most of which are cirrus passing by as verified by wide-field 
cameras on the top of PLATO when the Moon is up. Overall there is a greater fraction of 
photometric conditions at Dome A compared to Mauna Kea. However, these conclusions 
are tentative due to the limited sky field of CSTAR, different photometric systems, and the 
fact that the Dome A data are from only 5 months during 2008. 



3.6. Absolute Flux from the Sun 



The Sun and the Moon have o f course a strong effect on the sky background, especially 
when they are above the horizon. iLiu et al.l (120031 ) modeled the relationship between the 
sky brightness and the phase and elevation of the Moon. The way in which the Sun and 
the Moon affect the sky background should be formally analogous except for the additional 
complication that the Moon exhibits phases. 



Referring to the model of ILiu et al.l ( 120031 ) in V band, we fit the sky background flux as 
a function of the solar elevation with the same formula for our i-band data, 



Sun 



a!0 be + c, 



where -Fsun is the sky flux (assuming only a solar contribution and a background), 9 is the 
elevation of the Sun, and a, b and c are constants. We fitted for the constants with a nonlinear 



http : //www . gemini . edu/ sciops/telescopes-and-sites/ observing-condition-constraints 
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least-squares method, using the images with good transparency and no contribution from 
the Moon. The result is 

F Sun = 4.906 x 10 6 x io 339e + 370 (2) 

and is plotted against the observations in Figure [6j We find that the sky is about as dark 
as it ever gets in z-band once the Sun is below about 13° from the horizon. 



3.7. Astronomical Twilight and Hours of Dark Time 

Astronomical twilight is defined as when the Sun is between 12° and 18° below the 
horizon. However, this definition is really only appropriate at low altitude observatories and 
in the V band. At the high elevation of Dome A, and in i band, the sky is as dark as it 
ever gets when the solar elevation reaches about —13° (see Figure E])- Table [T] summarizes 
the numbers of hours when the Sun is located below elevations of —18°, —13°, and 0°. If we 
require only that the Sun is below the horizon, Dome A has more "night time" that any of 
the other sites. This is because Earth's elliptical orbit around the Sun results in the total 
hours of night time in a year increasing with distance from the north pole. However, for 
"dark time" to occur, we require the Sun to be at least a certain elevation below the horizon 
(conventionally, —18°). There is less dark time at Dome A compared to temperate-latitude 
observatories, although the —13° solar elevation limit makes this difference less dramatic 
than it would be at —18°. While the total dark time is less, the long periods of continuous 
darkness are an advantage for many observations. Over 83% images were obtained during 
2008 when the Sun was 13° below the horizon. The discussion in the rest of this paper is 
based on these images. 



3.8. Absolute Flux from the Moon 



The model for the sky surface brightness caused by moonlight is more complex. In 
our analysis, we ignore the variati on in the l unar brightness caused by the changes in the 
Moon— Earth distance. Following iLiu et al.l ( 120031 ). the apparent magnitude of the Moon 
can be approximated by an empirical formula: 



V(R, $) = 0.23 + 51ogP - 2.51ogP($), 



where R is the Moon— Earth distance, and $ is the Moon phase angle, and -P($) is the 
function of the full Moon illuminance. Then, using the same approach as for the Sun, the 
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sky surface flux contribution by the Moon, -Fmoon, can be expressed in the form of Equation 
flU multiplied by P($). Then, 

Fmoon = AP($)10 se + C, (3) 

where G is the elevation of the Moon, and A, B, and C are constants. C is the intrinsic 
background flux to be fitted. We applied these equations to the i- band data. Although the 
model is somewhat crude, it works well in the estimation of sky brightness illuminated by 
the Moon. For a m ore refined but slightly more c omplicated sky brightness model, we can 



refer to the study of iKrisciunas fc Schaeferl (Il99ll ) 



Figure [7] displays the Moon contribution to sky surface brightness as the function of the 
Moon phase and elevation. The sky brightness in the figure is the Sun-corrected sky surface 
flux without any contribution from clouds. The colored curved surface is the fitted model 
described by 

Fmoon = 300.576P($)10 0162e + 306.813. (4) 

We can see that the sky was brightest when the Moon was almost full. The outliers should 
be the points when aurorae were present. 

We compare the modified sky background (removing both the solar and lunar contribu- 
tions) with the raw one (before removing the solar contribution) in Figure [SJ The modified 
sky background becomes flat wherever the Moon is up. Figure [9] compares the histograms 
of the raw sky brightness in unit of mag arcsec -2 (the left panel shown in Figure [2]) with 
the modified sky background. The histogram with black slashes represents the original data 
and the one patched with red backslashes is for the modified sky brightness. The median 
sky brightness becomes 20.04 mag arcsec -2 after being corrected. About 80% of the images 
are fainter than 19.35 mag arcsec -2 . 



3.9. Correlation Between Sky Brightness and transparency variation 

The sky will be brighter when there is scattering of the Sun and the Moon light by 
clouds. In order to study this effect, we choose a day, June 16, when the Moon was nearly 
full, the Sun was far below the horizon, and the transparency was varying. Figure [10] shows 
the variations of the sky brightness and relative transparency in this day. We can see that 
they are strongly correlative, and the plot of the sky brightness against transparency in 
Figure [TT] confirms that the relationship is linear (again, outliers are the data polluted by 
aurorae). The line of best fit and the derived linear equation are also shown in the figure. 

The sky brightness is proportional to both the extinction (here, relative transparency) 
and the total illuminance of the Sun and the Moon. The larger the extinction, the brighter 
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the sky background for a given total illuminance. The greater the total illuminance, the 
brighter the sky background for a given extinction. 

We construct a model describing the sky brightness as a function of the total illuminance 
and transparency variation. When the relative transparency is negligible, the sky background 
should be the intrinsic sky background (a constant value), resulting in the model not con- 
taining an individual term of the Sun and the Moon. At new Moon or when the Moon is low 
enough below the horizon and the elevation of the Sun is low, the sky brightness should be 
the sum of the intrinsic background and cloud contribution. This requires a single relative 
transparency term. Then: 

-^sky = a (-^sun + E moon )E + bE + c, 

where F^y is the corrected flux of the sky light, F snn and -F moon are the model fluxes of 
the Sun and the Moon given by Equations (j2j) and (jlj) (without constant terms), E is the 
relative transparency and, a, b and c are constants, c denotes the fitted intrinsic background 
and b should be less than zero, indicating that the sky is dimmed by the clouds when it is a 
moonless night and the Sun is very low. The result, which gives an excellent fit to the data, 
is 

F gky = 5.4751(F sun + F moon )£ - 93.98E + 291.96. (5) 



3.10. Percentage of Images Affected by Aurorae 



For our SDSS i band (effective wavelength: 780nm; bandwidth: 160nm (IZhou et al. 



2010af )). the emissions of aurorae within this band mainly contain the lines from molecular 
nitrogen and molecular oxygen in the lower atmosphere and those from atomic oxygen in the 
higher atmosphere. The typical green— y ellow auroral emissio n lines are dominated by the 



atomic oxygen, lying outside the i band (IDempsey et al.l 120051 ) . To estimate the percentage 



of images affected by aurorae, we use the model of sky brightness as described in Equation 
flSJ). Figure [T21 shows how well the model fits the data and identifies the outlying points 
(implying the presence of aurorae) as circles. We define an outlier — and hence the image 
affected by an aurora — as a point that is more than 3a brighter than the model would predict 
(although it is also possible that a small fraction of those images are affected by airglow). 

There are about 3800 images polluted by aurorae, or less than 2% of the total. Most 
of the detected aurorae should be comparatively strong ones. The aurorae with weaker 
intensities are immersed in the fluctuation of the real observed data and errors of the fitted 
model. Therefore, the fraction of aurorae is a lower limit. However, we believe that any 
aurorae possibly present in the remaining 98% of our images are sufficiently faint as to not 
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affect photometric accuracy. We do not find obvious dependence of the aurorae on the 
transparency variation or cloud coverage, which implies that the above aurorae are detected 
reasonably, because the transparency variation distribution of those polluted images are 
generally uniform as we check. Figure [HI] shows the time distribution of the images affected 
by aurorae. A considerable number of them concentrate before or after 7:00 in UT when 
the Sun is highest. We can also find that there are few or no images affected by aurorae at 
about 20:00 nearly when the Sun elevation is lowest at 19:00. 



An optical spectrometer called Nigel ( iKenyon et al.l 120061 ) was installed in PLATO at 
Dome A in 2009 January and this will allow unambiguous identification of aurorae. 



4. CONCLUSION 

In 2008, CSTAR was deployed at Dome A in Antarctica to measure the night sky 
background, the weather conditions, and to study variable stars. About 310,000 images 
were acquired during the observing period of over 4 months. We have used these data to 
quantify the z-band sky brightness, variations in sky transparency, the background flux from 
the Sun and the Moon, and to derive some statistics on aurorae. 

Due to the high ecliptic latitude of the observed South Celestial Pole area, the sky 
background is little affected by the zodiacal light. The i-band sky brightness has a median 
value of 20.5 mag arcsec -2 when not affected by the Sun and the Moon. When corrected for 
the diffuse Galactic background, the sky brightness drops a further 0.06 mag. 

Comparing each image to our "best" , or reference, image observed on a clear and moon- 
less night, variations in transparency can be calculated. We find that the weather was stable 
and there was little cloud during most of the observing time. Conditions became somewhat 
worse and the transparency changed considerably in June and July when the Sun was at its 
lowest below the horizon in the year. Among the 83% images in which the elevation of the 
Sun was lower than —13°, about 67% of them were taken under conditions of little or no 
cloud, using the same criterion for the cloud cover as used by the Gemini Observatory at 
Mauna Kea. 

On the whole, Dome A has excellent optical sky backgrounds and had low cloud fractions 
during 2008. About 2% of our images were affected by relatively strong aurorae. 
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Fig. 1. — Sky brightness in i band (black dots). The red and blue curves are the elevations 
of the Sun and the Moon, respectively. 
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Fig. 2. — Left panel shows the histogram and the cumulative distribution function (CDF) 
of the z-band sky brightness distribution at Dome A during 2008. Right panel is the same 
information for the subset of images taken on moonless clear nights in 2008 June. 
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Table 1: Hours Per Year at Different Sites When the Sun is Below a Given Elevation. 



Site 


Latitude 


Elevation (m) 


Hours (< -18°) 


Hours (< -13°) 


Hours (< 0°) 


Dome A (Antarctica) 


-80°22' 


4093 


1680 


2606 


4471 


Mauna Kea (Hawaii) 


+ 19°50' 


4194 


3395 


3689 


4394 


Cerro Pachon (Chile) 


-30°14' 


2722 


3330 


3640 


4404 


La Palma (Spain) 


+28°46' 


2332 


3313 


3623 


4381 



Date(2008) 

Ma r. 14 Apr. 3 Apr. 23 May. 13 Jun. 2 Jun. 22 Jul. 12 Aug . 1 




I L_l i i i i i i I 

5.454 5.456 5.458 5.46 5.462 5.464 5.466 5.468 



JD- 2400000 x1Q 4 



Fig. 3. — i-Band sky transparency variation (in magnitudes) vs. Julian Date (JD) for all 
the CSTAR images during 2008. 
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Fig. 4. — (a) i-Band sky transparency variations (in magnitudes) from April 15 to April 25. 
The gaps are periods when sunlight made the sky too bright to observe, (b) May 10 to May 
20, (c) June 5 to June 15, and (d) July 15 to July 25. 



Table 2: Cloud Cover Conditions at Dome A. 



Flux 


Excess Extinction (mag) 


Fraction 


Cloud Cover 


< 50% 


> 0.75 


9% 


Thick 


50%-75% 


0.31-0.75 


17% 


Intermediate 


75%-90% 


0.11-0.31 


23% 


Thin 


> 90% 


< 0.11 


51% 


Little or none 



Note. — The first column is the observed flux of our standard stars relative to the reference image. The 
second column gives the z-band magnitude change (excess extinction) corresponding to the "flux" . The last 
column gives the description of the cloud cover. 
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Fig. 5. — Histogram of the number of images with the given transparency variation (in 
magnitudes) relative to the reference image. The curved line marked '+' is the CDF. 
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Fig. 6. — Relation between the sky brightness and the elevation of the Sun. The data points 
are those images that were taken in good conditions with negligible contribution from the 
Moon. The curved line is the fit from Equation (T5]). The dashed line at —18° gives the 
dividing elevation where the astronomical darkness is defined to begin. The dashed line at 
— 13° indicates a more suitable cutoff for i-band dark time at Dome A. 
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Fig. 7. — Relationship among the sky brightness (with the solar contribution subtracted and 
under conditions of good relative transparency), the lunar phase, and the elevation of the 
Moon. The curved and colored surface is fitted by the model in the form of Equation ([3]). 
The color bar is related to the Moon phase in percentage, which is the illuminated fraction 
of the full Moon disk. The outliers may be aurorae. 



-23 - 



Date(20QB) 




5.456 5.458 5.46 5.462 5.464 5.466 5.468 

JD - 2400000 4 

x 10 



Fig. 8. — Sky brightness before and after correction (subtracting the Sun and Moon contri- 
butions) and the contributions of the Sun and the Moon in the fitted models. The model 
of the Moon is displayed in continuous JD, but the model of the Sun is plotted at the time 
samples those same as the sky brightness. The original sky brightness and the model fluxes 
of the Sun and the Moon are shifted vertically from their proper values. 



Table 3: The Comparison of Cloud Cover Between Mauna Kea and Dome A. 



Mauna Kea (Gemini) 


Dome A 


Cloud Cover 


Extinction (V) 


Fraction 


Fraction 


Any other usable 


> 3 


10% 





Cloudy 


2-3 


20% 


2% 


Patchy cloud 


0.3-2 


20% 


31% 


Photometric 


< 0.3 


50% 


67% 



Note. — The definition of cloud cover is adopted from the Gemini Observatory. For comparison, we use 
V — i — 0.07 in extinction for the different tranparencies of these two bands as presented in the text. Note 
that the term 'photometric' as used here is just one kind of cloud cover category and it is different from the 
normal term 'photometric night'. 
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Fig. 9. — Histograms of the sky brightness distributions of both the original data and 
the modified one (subtracting the Sun and Moon contributions). The area patched with 
black slashes and the curved line with markers '+' represent the raw sky brightness and its 
CDF. The area with red backslashes and the curved line with circles are the modified sky 
background and its CDF. 
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Time (hour) 

Fig. 10. — Variations of the sky brightness and transparency variation (excess extinction in 
magnitudes) for 1 day in 2008. The black asterisks are the excess extinction and the red 
circles are the sky background. The x-axis is the time in hours elapsing from 2008 June 16 
00:00 UTC. 
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Fig. 11. — Correlation of the corrected sky brightness and transparency variation (in magni- 
tudes). The data points are observations taken on 2008 June 16, when the average elevations 
of the Sun and the Moon are — 23.°17 and 25.°12 respectively, and the lunar phase was ~ 95% 
($ = 25.8, -P($) = 0.54). The straight line is the fitted relation with the equation shown. 
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Fig. 12. — Observed sky brightness corrected for the Sun and Moon contributions vs. the 
values predicted by the fitted model from Equation flH]). The diagonal denotes where the 
observations and model agree exactly. The two dashed lines are 3a above and below this level. 
The circles show the observed data points, presumably affected by aurorae, with brightness 
brighter than 3a above the model. The transparency variations are in magnitudes. 
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Fig. 13. — Time distribution of the images taken during 2008 that were affected by aurorae. 
The counts are corrected for the image-sampling time distribution. 



